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Abstract 

<N ■ 

The noncommutative selfdual model in 6 dimensions is quantized and es- 
sentially solved, by mapping it to the Kontsevich model. The model is shown 



to be renormalizable and asymptotically free, and solvable genus by genus. It 
requires both wavefunction and coupling constant renormalization. The exact 
("all-order") renormalization of the bare parameters is determined explicitly, 
which turns out to depend on the genus sector only. The running coupling 
constant is also computed exactly, which decreases more rapidly than pre- 
dicted by the one-loop beta function. A phase transition to an unstable phase 
D ! is found. 

43 



1 Introduction 

This paper is the third part of a series of papers [1,2] studying the quantization of 
a noncommutative "self-dual" 3 model, by mapping it to the Kontsevich model. 
The model is characterized by an additional potential term in the action, which 
takes care of the UV/IR mixing following [3-6]. In the previous papers we discussed 
the cases of 2 and 4 dimensions, and showed that the model is renormalizable and 
essentially solvable. These cases are in some sense simpler because they are super- 
renormalizable. This is no longer the case in 6 dimensions, where new renormaliza- 
tion is required at each order, and the full complexity of an interacting quantum 
field theory is found with both wavefunction and coupling constant renormalization. 
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Indeed, the commutative 3 model in 6 dimensions is known [7] to be asymptotically 
free. 

Generalizing [1,2], we show in this paper that the self dual NC 3 model in 6 
dimensions can be renormalized and essentially solved in terms of a genus expansion. 
This is possible again using the results of [8, 9] on the Kontsevich model, which 
most however be properly renormalized. The model has 6 relevant parameters. In 
particular, both coupling constant and wavefunction renormalization are required in 
the 6-dimensional case, in addition to tadpole and mass renormalization which were 
sufficient in 2 and 4 dimensions. Remarkably, we find again that the renormalization 
is determined by the genus sector only, as in the case of 2 and 4 dimensions. 

After renormalization, all n-point functions can be computed in principle in terms 
of a genus expansion, and we give explicit expressions for the 1-, 2-, and some 3-point 
functions. We also find a critical surface a = in moduli space, which separates 
2 different phases. This provides a model which contains the full complexity of 
renormalization of a non-super-renormalizable asymptotically free quantum field 
theory, while being solvable and hence fully under control. It is very remarkable 
that a 6-dimensional interacting NC field theory allows such a detailed analytical 
description. It can therefore also serve as a testing ground for various approximation 
methods, which is of interest also in a more general context. 

In particular, we are able to determine exactly the RG flow of the bare parame- 
ters, as well as the running of the "physical" coupling constant i.e. the 1PI 3-point 
function. We find that the 1-loop beta function for the coupling constant correctly 
predicts asymptotic freedom, but wrongly predicts a (logiV)" 1 / 2 dependence on the 
scale as opposed to the correct (logiV) -2 dependence. 

Beyond the selfdual case Q — 1, perturbative renormalizability of similar models 
has been established using variants of a renormalization group approach, see e.g. [3- 
5,10-12]. The present approach using matrix model methods is certainly appropriate 
for the case Q — 1. We indicate in section [5] how it might be extended to Q ^ 1, 
which however is not carried out here. For a related matrix approach to NC field 
theory see also [13,14]. 

Being a continuation of our previous work on the 2- and 4-dimensional case, we 
will be brief in certain issues which have already been discussed there, and which 
apply without change. Nevertheless, the present paper is essentially self-contained. 
In section [5] we define the 4> 3 model under consideration, and rewrite it as Kontsevich 
model. We then briefly recall the most important facts about the Kontsevich model 
in section [31 The main technical analysis is contained in section 14.11 while the main 
results of this paper are collected in section H~2l The asymptotic behavior of the 1-, 
2- and 3-point functions is determined in section 14.31 and further aspects such as 
relation with string field theory are briefly discussed in section The perturbative 
results such as the one-loop beta-function are given in section El and we conclude 
with a general discussion and outlook. 
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2 The noncommutative selfdual (f> model 

We consider the noncommutative 3 model on the 6-dimensional quantum plane Mg, 
which is generated by self-adjoint operators^] satisfying the canonical commutation 
relations 

[xj, Xj\ iOij, (1) 
for i, j = 1, ...,6. We also introduce 

assuming that is nondegenerate. The model to be studied is defined by the action 

s= f + ^-<p 2 + n 2 {x4){xi4>) + ^- ( 3 ) 

Jr6 11 3! 

An additional oscillator-type potential fl 2 (x i (f))(x i (j)) is included following [3-6], mak- 
ing the model covariant under Langmann-Szabo duality, and taking care of the 
UV/IR mixing. The dynamical object is the scalar field = which is a self- 
adjoint operator acting on the representation space 7i of the algebra (CQ). The action 
is written with imaginary coupling iAo, so that the quantization is well-defined for 
real A ; otherwise the action would be unbounded. We will see however that after 
quantization, one can perform an analytic continuation to real Aq = 

Noting that the <9j are inner derivatives dif = — z[xj, /], the action can be written 

as 

/2 
-(xi(f)Xi(j) - XiXi4>4>) + Vt 2 Xi(pXi(t) + —<p 2 + 3 (4) 

using the cyclic property of the integral. For the "self-dual" point Q — 1, this 
simplifies further to 

S = J (f A + ^)0 2 + ^ 3 • (5) 

In order to quantize the theory, we need to include a linear (tadpole) counterterm 
—Tr(i\)A(j) to the action (the explicit factor iX is inserted for convenience). Re- 
placing the integral by J — > (2ir6) 3 Tr and adding a constant term for convenience, 
we are led to consider the action 

S = Tr^-Jcf + ^ 3 - (*A)A0 - ^J 3 - J A). (6) 

Here J = 2Z(27r6) 3 (^2 i x i x i + ^-) is essentially the Hamiltonian of a 3-dimensional 
quantum mechanical harmonic oscillator. A wavefunction renormalization Z has 

2 We ignore operator-technical subtleties, since the model will be regularized using a cutoff N. 



3 



also been introduced, which for the cubic term is absorbed in the redefined coupling 
constant A, in order to simplify the notation. The field will be the renormalized, 
finite physical field. In the usual basis of eigenstates, J then diagonalizes as 

J\ni,n 2 ,n 3 ) = 1Qtt 3 9 2 Z (m + n 2 + n 3 H )|m, n 2 , n 3 ), ^€{0,1,2,...} 

(7) 

assuming that 9ij has the canonical form 9 12 = — #21 — #34 — — #43 = ...=: 9. To 
simplify the notation, we will use the convention 

n = (ni, n 2 ,n 3 ), n = ni + n 2 + n 3 (8) 

throughout this paper, keeping in mind that n denotes a triple-index. It turns out 
that we need 

A = a + a x J + a 2 J 2 . (9) 
By suitable shifts! in the field 6 

= + ^-J = X + ^M (10) 

one can now either eliminate the linear term or the quadratic term in the action, 

S = T r(-^-M 2 4>+ l 4^) = Tr(-MX 2 + '4 X 3 - ^—M 3 ) (11) 
V 2i\ v 3! v J \2 3! 3(i\) 2 J v ; 

where 



M = ^.P + 2(i\) 2 A = xZVJ 2 + d (12) 

J = (13) 
Z x A Z 

x = v/l + 2(iA) 2 a 2 , (14) 

- - -(^r+^. w 

(fTTj) has precisely the form of the Kontsevich model [9]. The linear coupling of 
the field to the source M 2 resp. J 2 will be very useful for computing correlation 
functions. J now has eigenvalues 

Q -4- 3 

J\n 1 ,n 2 ,n 3 ) = 16tt 3 9 2 (m +n 2 + n 3 + — ^— )|ni,n 2 ,n 3 ) (16) 

which will be finite after renormalization, and 

^ = ^-^ = -15^^- (17) 



3 for the quantization, the integral for the diagonal elements is then defined via analytical con- 
tinuation, and the off-diagonal elements remain hermitian since J is diagonal. 
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2.1 Quantization, partition function and correlators 

The quantization of the model (EJ) resp. (TTTT) is defined by an integral over all 
Hermitian iV 3 x N 3 matrices 0, where N serves as a UV and IR cutoff. The partition 
function 

Z(M) = J D0 exp ( - Tr( - ^M 2 4> + ^ 3 )) = e F ^ (18) 

defines the "free energy" F(M), which is a function of the eigenvalues of M resp. 
J. Since N is finite, we can freely switch between the various parametrizations ([6]), 
ffTTj) involving M, J, 0, or 0. Correlators or "n-point functions" are defined through 

{<t>ixh-4i n 3 n ) = tt(m) J D ^ ex P(- S )^ii-4j„, (19) 

keeping in mind that each index denotes a multi- index ([8]). Using the symmetry 
Z(M) = Z{U~ l MU) for U e U{N 3 ), we can assume that M is diagonalized with 
ordered eigenvalues m^. There is a residual U(l) x U(3) x U(6) x ... symmetry, 
reflecting the degeneracy of J. This implies certain obvious "index conservation 
laws", such as (<f) kl ) = 5 kl (4> u ) etc. 

The nontrivial task is to show that all these n-point functions with finite indices 
have a well-defined and nontrivial limit N — > oo, for a suitable scaling of the bare 
parameters. In addition, the index dependence of these n-point functions must 
be nontrivial. Recall that the wavefunction renormalization is already taken into 
account in ([7]), so that <fi is the finite, physical field. The free parameters should in 
principle be determined by choosing renormalization conditions, such as (0oo0oo) = 
5tt ^ fl+i ' e ^ c ' These conditions can easily be solved using the explicit results given 
in sections 14.11 and 14.21 relating the bare parameters with the "physical" quantities. 

Noting that the field couples linearly to M 2 resp. J 2 in ffTSl) . one can compute 
the connected n-point functions by acting with the derivative operator 2i\-^ — 
2«^^2^jj on F. Anticipating some results below, we introduce the quantities 

iX * = Jz~^ r = ^z (20) 

which turn out to be finite after renormalization. Then 

(kn-^jJc = {2i\ R - w -)....{2i\ R - w -) F(J 2 ). (21) 

Since the connected n-point functions are independent of the shifts = 0+^ (|TU|) 
for n > 2, the lhs coincides with the desired correlator (<p il j 1 ...(f) in j n ) c for n > 2. This 
strongly suggests that \r should be finite. 
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Using the Dyson- Schwinger equations for the path integral (JT51) . one can derive a 
number of nontrivial identities for the n-point functions. Since their derivation was 
already given in [1], we simply quote them here with the appropriate substitutions. 
In particular, one finds for the propagator 

(4>ki(l>ik) = o (<Pkk - <Pu) (22) 

m k ~ m i 

for k 7^ I (no sum), where denotes the eigenvalues of M. This gives 

{<pki<Pik) = - 2r - + ^ XR j 2 {<Pkk - <t>u) (23) 
Jk + Ji Jk~ J i 

using (112p . The first term has the form of the free contribution, and the second is 
a quantum correction. The latter reduces to the 1-point functions, which will be 
obtained from the Kontsevich model. Similarly, one can show [1] 



2i\ 



>kl9lk(Pkk } = —r, V \(Pkk - (Pll 9kk - 9kl9lk 

m k ~~ m l 



(24) 



(no sum) for k ^ I. Using (l2"3"j) . this gives the connected part 

2iX 



ki4>ik4>kk) c = — o— — 2 ((((fikk - 4>u)4>kk)c - (4 > ki4>ik)') ■ (25) 
m k ~ m t ^ ' 



Clearly these relations can be generalized, greatly reducing the number of indepen- 
dent correlators. However, we will establish finiteness of the general correlation 
functions directly, by showing that the appropriate derivatives of the generating 
function F(J) are finite and well-defined after renormalization. 



3 The Kontsevich model: facts and background 

The Kontsevich model is defined as a matrix integral 



?Kont( »>\ F Kont 



f dX exp j 




'} 


f dX exp j 


[-Tr(Sf) 


} 



(M) = e F = "—^ , " /J (26) 



over hermitian iV 3 x A^ 3 matrices X, where the parameter M is some given her- 
mitian iV 3 x A^ 3 matrix. This model has been introduced by Kontsevich [9] as a 
combinatorial way of computing certain topological quantities (intersection num- 
bers) on moduli spaces of Riemann surfaces with punctures. It turns out to have 
an extremely rich structure related to integrable models (KdV flows) and Virasoro 
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constraints. For our purpose, the most important results are those of [8,9,15] which 
provide explicit expressions for the genus expansion of the free energy. Note that A 
can be introduced by rescaling the variables, writing M = X 2 ^ 3 M. 

The matrix integral in ( 1261) and its large N limit can be defined rigorously in 
terms of its asymptotic series. This involves the variables [9] 

t r : = -(2r + 1)!! 02r+i> Or ■= \TrM- r . (27) 

Then the large N limit of the partition function Z Kont (M) can be rigorously defined 
[8,9], which turns out to be a function of these new variables only, Z Kont (M) = 
Z Kont (8i). In our case, 8 r is divergent for r = 1, 2, 3, which indicates that the model 
requires renormalization. 

One can furthermore consider the genus expansion 

^Kont pKont y pKont (28) 

9>0 

by drawing the Feynman diagrams on the appropriate Riemann surface, as usual for 
matrix models. In [8] it was shown that all F^ ont can be computed using the KdV 
equations and the Virasoro constraints, which allows to find closed expressions for 
any genus g. They are given in terms of the following variables 



hfaU) = £ = -(2fc-l)"E ,- 2 I 7^T > ( 29 ) 



where u is given by the solution of the implicit equation 

wo = Iq(uo,U) = -V 7=== (30) 
i V m i - 2u o 

These variables turn out to be more useful for our purpose. Using the KdV equa- 
tions, [8] obtain the following explicit formulas: 

Ft nt = JE^^E^- 2m o) 3/2 - u oE^- 2w o) 1/2 

i i i 

I U ° 1 V In f ( ™' ? - 2Uo)1/2 + {Tnl - 2Uo)1/2 ) (31) 
6 2^ 1 rhi + fhk ' 



i . fc 



F, Kont = — ln^^, (32) 
1 24 v ; 

«r* = iL [ 5 <r^ + 29 n^h + 28 (dy • < 33 > 
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etc. This form of F was first found in [15]. The sums over multi-indices (jHj) are to 
be interpreted as 

N-l 
i ii,i2,i3=0 

truncating the harmonic oscillators in (J^J). For our purpose, the most important re- 
sult is that all F^ ont with g > 2 are given by polynomials F^ ont = ^2 X{i k } xl £--- x ig-2 
in Xk = ^jfc/(l — h)^ 1 with the constraint 

(k-l)l k = 3g-3, (34) 

2<fc<3g-2 

where is the power of x^, and XlM is some (rational) intersection number. 

While many of these expressions are divergent as iV — ■> oo for J given by f fTST) . 
the physically relevant observables will be convergent after renormalization. 



4 The Kontsevich model applied to the 3 model 

In order to apply the above results to the noncommutative 3 model, we need the 
following slightly modified version of the Kontsevich model, corresponding to the 
action (ITT]) : 

Z(M) = exp(F(M)) = Z Kmt [M]Z free [M] exp(— TrM 3 ) (35) 

o 1 1 A ) 

where 

F := F Kont + F free + ^l^TrM 3 

i i 

and F fl = F^ ont for $ > 1. In the present case, the eigenvalues rhi are given by (FEZ]) 

rhi = X~ 2/z xZ\J jf + d. (37) 

The model as it stands is ill-defined for iV — > oo, since m and many of the above 
sums are divergent. However, we can recast (I36p using more appropriate variables, 
which suggests how to renormalize the various parameters, i.e. how to scale them 
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with N. Note that only the combinations a/ mf — 2w enter in (I2"9"j) and (1361) . which 
can be rewritten using (fT2l) as 



m 2 - 2u = a;ZA~ 2/3 a/ J 2 + 26. (3* 



Here 6 is defined through 

x 3 Z 3 / 1 



fe--d)=xZA- 2/3 ii = -^ 1 (39) 



using the constraint (1301) . which is replaced by (1391) henceforth. Eliminating wo using 
( 1391) and expressing A in terms of Xr, F takes the form 

1 HTr, T 3 b-U 



3(zAfl) 2 a;Z 4^ v («Aij) 2 xZ 4^ V 1 



(6-±d) 3 1_ 1 



Q(iX R ) 4 x 2 Z 2 2 _ 



The quantities J, b and Xr will be finite after renormalization, rendering the model 
well-defined. We consider F = F(J 2 ) as a function of (the eigenvalues of) J 2 from 
now on, while b is implicitly determined by f[3"9"l) . Since the eigenvalues </& only enter 



through the combination y j| + 26, we note that the eigenvalues can be analytically 
continued as long as this square-root is well-defined. 

We can now compute various n-point functions by taking partial derivatives of 
F = J2 g Fg with respect to J 2 , as indicated in section 12.11 For the "diagonal" n- 

point functions (<j>u--.<j>kk)c, this amounts to varying the eigenvalues Jf. In doing 
so, we must remember that b depends implicitly on J 2 through the constraint (}3"9"1) . 
However, some of these computations simplify recalling that the constraint (1391) for 
b arises automatically through the e.o.m: using 



Jf + 26 



we can write 



%-MJi) = ^-MJl b) + JU(Jf ! 6)4-6 = -2-F (J 2 ; 6). (42) 



dJ 2 uv w dJ 2 uv 11 ' db uv l ' 'dJ 2 dJ 2 

Thus for derivatives of order < 2 w.r.t. J%, we can simply ignore 6 and treat 
it as independent variable, since the omitted terms (I4ip vanish anyway once the 
constraint is imposed. 
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4.1 Renormalization 



This section contains some computations required to show finiteness of the n-point 
functions for any genus g. This will establish Theorem 1 in section 14.21 

The 1-point function We can now determine the required renormalization by 
considering the one-point function. Using (HQ]), (I4B and (1391) . the genus zero con- 
tribution is 



■>kk) g =o = 2zA flT ^ F (J 2 ) - 



9 v (p\_ Jk 
dJl W } i\ R x*Z* 



1 , ^ 1 r 



^yk + (i^R)y~\ , — J k + i\ R Z 2 ai 



.1 



i y k J Jf + 26 + ( J 2 + 26) <A « 



=: W(j/ fc ) (43) 
using ( l20i) . which must be finite and well-defined as N — > oo. Here we define 



y k = pl + 2b. (44) 

We will find that W(y) becomes a smooth function of y after renormalization, which 
amounts to the statement that the index dependence of the one-point function is 
"smooth". This is a typical feature of matrix models, reflecting some "smooth" 
distribution of eigenvalues. This becomes here part of the statement of renormaliz- 
ability. To proceed, we need to understand the function 

f(y) ■= £ / 1 ( 45 ) 

3 2/^ + 26+ (J 2 + 26) 
which as it stands is ill-defined for iV — > oo. 
4.1.1 Renormalization of f(y) 

In order to make sense of f{y), we consider the Taylor-expansion of 

f(y; J) = /Z i 1 = MJ) + vW) + Sr(v\ J) ( 4 6) 

i 2 / v /j| + 26+(J| + 26) 
in y, in analogy to the usual strategy in renormalization. Here 
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f R (y;J) = y 2 £ ; -. (49) 

i (Jj + 2b) (yyjjj + 2b + {J] + 2b)] 



where /i, f 2 are divergent constants obtained by fixing J as in (TTBT) . while the fi^SJ) 
are regular (convergent) functions obtained by taking into account variations SJ of 
the eigenvalues. This is necessary e.g. to compute partial derivatives w.r.t. J. Thus 

1 f N 1 

f° = fia 3zm2 / dx i dx 2dx 3 =7-5 + finite 

(167T^2)2 ^ ^ + ^ + X3 + 3+^ )2 

= ^^((61og(2)-31og(3))iV-^Mlog(iV))+Mzte (50) 
and similarly 

h = 4(i6^W log(iV) +/zmte - (51) 

The remaining part ,/«(?/; J) is well-defined and convergent, provided b and i.e. 
/iij are finite, which will be assumed from now orfl To understand it better, we note 
that 

f'M = f"(y) = 2j2 J l + 2b 



(y^jl + 2b+(J* + 2b) 

is positive, and similarly f' R (y) > 0. Hence fn(y) is a rather simple convex smooth 
function of y > 0, which satisfies /k(0) = f' R (0) = 0, and it remains only to 
determine its asymptotic behavior for large y ~ Jjy, i.e. small x. This can be 
obtained by writing 

v r^N/y i 

f R (y) « n JLm 3 / rf3x 7 ; v ( 52 ) 

( 167r r iJo/j/ (x 2 + 26/y 2 ) ( v 7 ^ 2 + 2&/V + + 26/y 2 )) 

where x = ^. This integral is convergent for large x, but logarithmically divergent 
at the origin i.e. for large y ~ J^. We can hence replace the upper integration limit 
in (152j) by x = 1, to obtain the asymptotic behavior for large y: 

y f i ..11 y 



/*(*) ^ (16^ J Jo/y dx ( x i + 2b/y^ W 2 (16^2)3 J,1 ° g( ^ ) - (53) 



'We do not pursue the possibility of divergent which does not appear to be interesting. 
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4.1.2 Renormalization of {<pkk)g=o 

We have seen that only f and fx are divergent in (H6l) . while f R (yk) is finite and 
well-defined as N — ► oo. Then fl43l) becomes 



Ws=o = tt^T^ ~ 77TT + ^A i? Z 2 ai + iX R (f (J) + Vk fi(J) + /r(^; ■/)) 



iX R xZ 

\ R xZ ' J aK i\ R 

+iX R (f R (y k ; J) + fo, R (SJ) + VkfxA^J) 



(i\ R Z 2 a x + iAfl/o) + f 1 + iAjj/i) y fc 

\lAnXZi / 



c _ _ rJ, + + / Q + y k f ltR (5J)) (54) 



Here we define the quantities 

xZ = 1 

c = iX R (Z 2 ai + f ). (56) 

Since Jk, yk and f R (yk) are independent functions of (as long as b ^ 0), it follows 
that (4>kk) g =o is finite if and only if the four quantities 

1 r 

(iX R , c, — — , — ) (57) 
i\ R a i\ R 

are finite. Remarkably, the condition that a be finite will also guarantee that the 
higher genus contributions are finite. The second form of (I55I) implies that this 
is possible only for real coupling (iX R ), since x and Z should be positive, while 
fx ~ -logiV. 



4.1.3 Higher derivatives and higher genus contributions 

The connected part of the n-point functions for diagonal entries (4>i 1 i 1 ---4>i n i„) c are 
obtained by taking higher derivatives {iX R ) n -^...-^j r °f F{J 2 ) resp. ( 1541) . Since 

*1 in 

the (infinite) shift = 0+^ drops out from the connected n-point function for 
n > 2, these coincide with (<f>i 1 i 1 ---<j>i n i n ) c for n > 2, and we have to show that they 
are finite. 

Consider first the genus contributions. To compute higher derivatives of Fq 
w.r.t. J 2 , we must also take into account the implicit dependence of b on J 2 . Indeed 
b is a smooth function of J as shown in ( 166]) . Furthermore, recall that Mr Fo(J 2 ; b) 
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vanishes through the constraint (I4ip . however this is no longer true for the higher 
derivatives. In particular, taking derivatives of (14 ip we find 

92 rn , r2 n d f b-\d ^ 1 \ 1 - I x 



dfidb^*" a/3V(zA*)%Z 



+ 

which are both finite and smooth using f l62l) below. Combining this with the explicit 
form ( |54|) and using the results of section 14.1.11 it follows that all higher derivatives 
of F (J 2 ) w.r.t. J\ are finite. 

For the higher genus contributions, we also need the derivatives of the quantities 



In particular, 



h(Jlb) = -^Zj2 {j2 ^ 2b)l =X 2 R xZf 1 (J) 



= XlxZih + f^SJ)) (60) 

where fi,R(5J) is a finite and smooth function (l4"g]l of Jf and b (which vanishes for 
5 J = 0). Hence for Ii to be finite we need 

-±—^f x ~\ g{N), (61) 

which using implies that all higher vanish for N —>■ 00 unless (1 — ii) ~ 0. 
Hence the higher-genus contributions can be nontrivial only if we carefully take a 
"double-scaling" limitEl and require J * Wi to be finite, while I\ — > 1. Using ( 1601) . 

(l-Zi) - 3~ 

we have 

T5^(l " h) = -7^5 - /i^J) (62) 
a r xZ a{iAR) z 



where a is defined in ( 1551) . It follows that 

h (2k - 1)!! 



rE 



2fc+l . — _ 1 

~ -(2A;-l)!!((^A i? ) 2 a)~ r (63) 



5 This is quite reminiscent of the double-scaling limit of matrix models in the context of 2D 
gravity 
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(the last form holds for 5 J = 0) is finite and nontrivial for k > 2, provided 

a = finite, (64) 

assuming that iX^ is finite. In particular, all derivatives of F g for g > 2 w.r.t. J| 
are manifestly finite, and thus all genus g > 2 contributions to the diagonal n-point 
functions are finite. Using ( 162]) and ( 132|) . this also holds for genus 1. 
Finally, from the constraint ( 1391) we derive 

^ 2 ^r*/#7i 



7 J 2 + 26 i Jjl + 26 d ^ 



which using ( lfi"2"]) gives 



9,11 1 



This is again finite, using the above assumptions. Hence 6 depends smoothly on 
the variations in J*, provided a is finite. If a = 00, then the constraint (1391) cannot 
be solved any more for 6 as a function of the eigenvalues J (and their variation), 
rendering the model non-renormalizable. This implies that the genus sector fully 
determines the required renormalization, as was found previously in 2 and 4 dimen- 
sions [1,2]. 



4.2 Main result and renormalization group flow 

We have now established all required formulas for the derivatives of F(J), and hence 
for the diagonal n-point functions. We also showed that all of these are finite and 
have a well-defined as N — > 00, provided a few renormalization conditions hold. Let 
us collect these conditions and use them to determine the required scaling of the 
bare parameters. 

Assume first that Xr = 0. Imposing this exactly (i.e. independent of N) implies 
a = —xZ and iX = iXna 2 = 0. This is the free case with x and Z finit^, as can be 
seen e.g. from ( 1541) and ( |83j) . 

Hence assume A# 7^ 0. We assume also that 6 7^ (we will see that this restriction 
is in fact not important), so that the 4 quantities in (I57p are finite. By taking 
products, this implies that r and — are also finite, hence cr^O. Finite a is in fact 

6 We shall not pursue here the possibility of other more subtle scaling limits. 
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also required for the higher-genus contributions to be finite. Using (I5"5"j) . we then 
obtain 

h = ~ \ = los(N) ~ \ (67) 



which implied 



iX = iX R x 2 Z 2 = ^ -3 ~ (iA*)" 3 log(Ny 2 (6* 

i | i (i\ R )2\og(N) y 



2 

a "T" 2 (16tt 3 2 ) 3 J 

for large iV. Note that in a perturbative approach i.e. if this is formally expanded in 
terms of Xr, this is divergent and requires renormalization at each order. Neverthe- 
less the closed form ( 168|) is rather simple, with leading behavior (iX R )~ 3 log(iV) -2 . 
Z is then determined through 

Z = r^- ~log(iV)- 2 (69) 
iX R 

since r = is finite, and similarly 

x = — L= ~ log(JV). (70) 



'rZ 

This gives 02 using (|14p . a x is determined from 

Z 2 ai = ^--f (71) 

whose leading term is linearly divergent in N, and gives the mass renormalization 

2 2(»A fl ) 2 1 y2 „ 2{iX R f f 

^ e = -^w- r Za ^m^- (72) 

using (TT7T) up to finite corrections. Its leading dependence on N is given by ( l50i) . 
Finally, d is determined by the constraint 



rf = 2fr + 2(zA fl ) 2 xZ V— (73) 

'jf + 2b 



which involves the further finite parameter b. Then a follows from ffToT) . which is 
quadratically divergent in N in agreement with the perturbative result (11211) . 

Finally, (|54|) shows that a = marks some singularity or phase transition, divid- 
ing the moduli space into two disconnected components with a ^ 0. The explicit 



7 note also that (|67|) implies that (i\n) should be real, since xZ should be positive. 
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form of the 2-point function (!HH|) below shows that the phase with a < is the 
"physically relevant" one, while {4>ij(pji) < for a > for small indices Since 
the matrices are supposed to be Hermitian, this signals an instability or condensa- 
tion for the low modes, i.e. some kind of tachyonic behavior. Observe also that e.g. 
(0oo) = can only be realized for a < 0. Nevertheless, the renormalized n-point 
functions remain well-defined for a > 0. This phase transition is only seen in the 
genus sector. 

Finally consider the case b = 0, so that y k = J k . Then finiteness of the one-point 
function requires only that the three quantities 

(iA*, c, ^- + ) (74) 
iX R a iX R 

are finite. However, the genus zero result for the 2-point function (15H|) shows that 
- must also be finite, and we are back to the previous analysis with finite 

Putting these results together and recalling the structure of the higher genus 
contributions F g stated below ( 1331) . we have established the following: 

Theorem 1 All derivatives of F g w.r.t. J| for g > (as well as all functions 
F g for g > 2) are finite and have a well-defined limit N — ► oo, provided the 6 
quantities (iX R , a, fi R , r, c, b) are finite and fixed, and 

This determines the scaling of the 6 bare parameters 

(i\(N), Z(N),fi 2 (N), a (N), ai (N), a 2 {N)) (75) 

as a function of N through fl§ED, f3D, fZ2J, dZSD resp. pS]>, <|7T|>, and (f7D|). This 
defines a renormalization group flow in the space of free parameters. If desired, we 
could now impose some specific renormalization conditions such as (4>oo) = etc. 
These would then determine the renormalized parameters iX R etc., which in turn 
would determine the bare ones. 

The most interesting "real" sector is given by iX R G M, a < 0, r > 0, [i R > —3 
and b G K such that J| + 2b > for all k. Then iX G M for large enough N, and Z 
and x are positive. The case a = is a singularity of the genus sector. 

Since the connected n-point functions are given by the derivatives of F — 
5] >0 F 9 w.r.t. J, this implies that all contributions in a genus expansion of the 
correlation functions for diagonal entries (<f)kk---<f ) ii) are finite and well-defined. The 
general non-diagonal correlation functions are discussed in section 14.2. 2| and also 
turn out to be finite for arbitrary genus g under the same conditions. Putting these 
results together we have established renormalizability of the model to all orders in 
a genus expansion, i.e. 

Theorem 2 All connected genus g contributions to any given n-point function 
{(f>i 1 j 1 ---(t ) injn)c ore finite and have a well-defined limit N — > oo for all g, under 
the above conditions. 
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Moreover, they can in principle be computed explicitly using the above formulas. 
This immediately extends to non-connected diagrams (<f>hj l ---4 > injn) ■ Furthermore, 
( j95l) shows that any contribution to F g has order at least ({i\R) 2 a) 2g ~ l . This implies 
(but is stronger than) perturbative renormalizability to all orders in Xr. This might 
not seem surprising in view of the results in [3-5], however note that the present 
model is more complicated than the </> 4 model where the beta function was found to 
be zero at one loop for Q = 1 [16]. 

It is worth pointing out that only the genus contribution requires renormaliza- 
tion, while all higher genus contributions are then automatically finite. This is very 
interesting because the genus contribution can be obtained by various techniques 
in more general models, see e.g. [17]. A related approach was studied in [13, 14] 
without the oscillator potential. 

It may appear surprising to find a well-defined 3 model for real coupling, where 
the action is not bounded from below. This is possible because we first quantize 
the model for imaginary coupling, where the Kontsevich model is well-defined. The 
full genus expansion is then available, and the limit iV — > oo is under control upon 
proper renormalization. This allows in a second step to define the model for real 
coupling i.e. real iX^_, through analytic continuation. Similar behavior is well-known 
in the context of pure matrix models, which typically allow extension to seemingly 
unstable potentials [18]; this is usually interpreted as suppression of tunneling from 
a local minimum. 

4.2.1 Exact renormalization and fine-tuning 

The scaling (I68p of the bare coupling is derived assuming that the one-point function 
(IB"4"|) and therefore its coefficients a and Xr are kept fixed and independent of N. 
This is the usual way to proceed in perturbative renormalization, and it would lead 
to new infinite renormalization at each order of perturbation theory, as can be seen 
by expanding (1681) formally in terms of Xr. 

On the other hand, we have the full renormalization in closed form available, 
which shows a simple leading scaling law for large N. It is then natural to ask what 
would happen if we scale the bare coupling in a simpler way, e.g. respecting only 
the leading scaling 

iX = c x log(iV)- 2 , Z = c z log(iV)- 2 , x = c x \og(N) (76) 

for some constant C\,Cz,c x rather than the exact form ( 168]) etc. It turns out that 
this is not sufficient, and the renormalization must respect the exact form of ( |68l) ff. 
rather than just their leading terms. 

To see this, we determine the quantities (Xr, a, r, ...) for the simpler scaling (1761) . 
and check whether they also converge as iV — > oo. This is obviously the case for 
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A/? ~ and r ~ -4—, while a is given by (JHID 

I = __L _ (iAR)Vl = log(JV) + i^^^j log^ (77) 

For this to have a well-defined limit A" — > oo, the above scaling ( 1761) is not sufficient, 

2 



but the constants of proportionality must satisfy = ( 167r 3g2)3 ; and moreover 
the scaling must be refined to give the sub-leading constant (1771) . 

This shows that even this relatively simple, solvable asymptotically free model 
has a rather severe "fine-tuning problem" , i.e. there is no obvious naturalness in the 
scaling of the bare parameters. 

4.2.2 General n-point functions 

Finally we show that all contributions in the genus expansion (and therefore per- 
turbative expansion) of any n-point functions of the form 

(0ilil--0Wn) (78) 

have a well-defined and finite limit as iV — > oo, provided the above renormalization 
conditions hold. The argument is the same as in [2], which we repeat here for 
convenience. 

Recall that the insertion of a factor can be obtained by acting with 2i\^-A T 

on Z(J 2 ) , resp. F g (J 2 ). Now any given correlation function of type (1781) involves 
only a finite set of indices . Thus taking derivatives w.r.t. Jfj amounts to 

considering matrices J of the form 



J 



diag(Ji, ...J fc ) + SJ kxk 








diag(J fc+ i, ...Jjv) 



(79) 



where k is chosen large enough such that all required variations are accommodated 
by the general hermitian k x k matrix 

Jkxk ■= (diag(Ji, ... J k ) + 5Jfcxfc) (80) 

in ( 1791) . while the higher eigenvalues J k +i, ... Jjv are fixed and given by ( fTBT) . Therefore 
we can restrict ourselves to this k x k matrix, which is independent of N . As was 
shown in section I4.1[ all F g are in the limit — ► oo smooth (in fact analytic) 
symmetric functions of the first k eigenvalues squared, hence of the eigenvalues of 
(Jkxk) 2 - Such a function can always be written as a smooth (analytic) function of 
some basis of symmetric polynomials in the J 2 , in particular 

F g {Jl...Jl) = f g (Tr(J 2 kxk ),...,Tr(J 2 k k xk )). (81) 
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This can be seen by approximating the analytic function F g (zi, ..,z k ) at the point 
Zi = Jf by a totally symmetric polynomial in the Zi, which correctly reproduces the 
partial derivatives up to some order n. According to a well-known theorem, that 
polynomial can be rewritten as polynomial in the elementary symmetric polynomi- 
als, or equivalently as a polynomial in the variables s n := z™, n = 1, 2, k. This 
amounts to the rhs of flHTT) . 

In the form flHTl) . it is obvious that all partial derivatives -|g- of F„ exist to any 

given order, and could be worked out in principle. This completes the proof that 
each genus g contribution to the general (connected) correlators (<fii 1 j 1 ----<fii n j n ) is 
finite and convergent as iV — > oo. 

4.3 Asymptotic behavior of the correlation functions 
4.3.1 (<f> kk ). 

Using (1531) . the behavior of (<pkk) for large k is dominated by f^yk), so that 

2(16^F^ l0g( X 



tt ^^9 n^3 ^ lo g(T) forfc^oo. (82) 



4.3.2 (Mik). 

We can use (1231) to obtain the genus contribution to the 2-point function (4>ki4>ik) 
for k I. Using (loTj) we obtain the exact expression 

/, , v 2o; _1 ^fniVk) - fR(Vl) / q o \ 

{<Pkl<Plk)g=0 = ; \-2(tX R ) 2 ( 83 ) 

Vk + yi yi- vi 

Note that the free case corresponds to i\r = and a = — 1. Then indeed also the 
higher genus contributions vanish. 

Consider the behavior for large k « I, which for large indices is dominated by 
the terms involving f R : 

ii i \ _ /- x ^fRivk) - fR{yi) 

{(PkifPik) ~ 2(2A K ) 2 

~ (i^r) 2 — -i— fniVk) 
Vkdy k 

for large k m I. It is worth pointing out that this 2-point function (and similar the 
3-point function (190]) below, etc.) is essentially determined by a function of a single 
variable y, which describes the dependence of the one-point function (4>kk) on the 
index k. This is characteristic for the genus sector, which is essentially determined 
by an eigenvalue distribution. 
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4.3.3 (4>u4>kk)- 

As a further example, consider the 2-point function {<j>u<j>kk) f° r k ^ I, which vanishes 
in the free case. To compute it from the effective action, we need in principle 

(<Pu<Pkk)c = (4>u4>kk) - (kk)(4>u) = (2iX R ) 2 - 1 -- 1 -(F + F 1 + ...)• (85) 

dJf dJ£ 

Even though this corresponds to a nonplanar diagram with external legs, it is ob- 
tained by taking derivatives of a closed genus ring diagram. Therefore we expect 
that only F will contribute, and indeed the derivatives of Fx contribute only to 
order X R . For k ^ I, the genus zero contribution is 

d d 

{(t>u4>kk)c = (2i\ R ) 2 F 

dJf <9J£ 

= HXn) 2 1 



J 2 + 2bJj 2 + 2b \Jj* + 2b+Jj? + 2b 



2 



W I ] (86) 

as k ~ I — > oo, using again (jHJ). For k = I, we need 

(<Pkk<Pkk) c = (2iX R ) 2 -^-F » ^rJR^ SJ), (87) 

d{J 2 ) 2 2 dJ 2 k 

since fFt(y) is the leading contribution for large k. Writing 



( \ 



dJ 2 ' ' ^ ' v ' dJl 



[ (J 2 + 2b) (y ky Jj* + 2b + (J 2 + 26)) ) 



we note that the first term (involving the sum) is dominant, while the second is only 
one term in a large sum. This then gives the same asymptotic behavior as the usual 
propagator, 

(2i\ R ) 2 d_ 

9y 2 k 

for large k, using On the other hand, it is very different from {4>u4>kk)c given 

by (I86p . This should not be a surprise, since the latter corresponds to nonplanar 
diagrams as in figure [3 



kk4>kk) c ~ — 7T— -K-^fn{yk] 5 J) ~ {<t>ki<t>ik) (89) 
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Note that in (IHSj) . a derivative w.r.t. Jf has essentially been replaced by a 
derivative w.r.t. y = y^. This is again characteristic for the genus sector. In fact, 
Fo can be obtained from the Dyson-Schwinger equation by making precisely this 
approximation. 



4.4 The 3-point function 

We can similarly use the exact expression (1251) for the 3-point function to determine 
the asymptotic behavior of the 3-point function at genus 0: 

(4>ki(fiik4>kk) c = 2zA R ~ 1 ~ {{{4>kk - 4>u)4>kk)c - (<f>ki<f>ik)) 

J l ~ J f 

vl-yi 2 \ d yl Vk-yf 
i d 



2 Vk dy k 

^R?\ h , L 3 ^°g(f) (90) 
4 (I6ir d 9 2 y y% J 



for large k ~ I, using flHbj) . (1841) and floBl . The same behavior is found for 

d 3 1 <9 2 

(MkMc = (2i\ R ) 3 —~—F « (2iX R ) 3 -—^-f R (y k ), (91) 

up to a factor 2 which reflects the exchange symmetry between the external "legs" 
4>ki and 0ifc for k = I. Combining (1901) with (1541 and ( 1891) . the IPI vertex behaves 
like 

(4>kl4>lk4>kk) ip/ = 7~7 T \ 7~7 T \ "7~7 7 T" (4>kl4>lk4>kk) c 

{<Pkl<Plk) c {<Pkl<Plk) c (<Pkk<Pkk) c 

« -2(iA K )^ 3 (167r 3 6' 2 ) 6 f L^ 2 ^_(J_Y (92) 

for large k ~ I. This result is exact, because the higher genus contributions de- 
cay more rapidly as discussed in the next section. In particular, this establishes 
asymptotic freedom. 

Several remarks are in order. First we note that this vertex function decays like 
( l0 g fc )2 rather than ^ log ^i/ 2 , which would be found by a 1-loop computation as shown 
in sectional Moreover, (I9"2"j) is in nice agreement with (1681) . demonstrating that the 
IPI vertex approaches the bare coupling — y for k — > N. Finally, note that the 
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effective coupling constant for this asymptotic domain is ■A- rather than X R . This 
is clearly a purely "quantum" effect, which is again due to higher-order correction 
to the one-loop RG result (11 13j) . 



4.5 Higher genus contributions 

We illustrate here the stronger decay behavior of the higher-genus n-point functions, 
in the example of the one-point function. Consider 



kk/g=l 



2i\ 







R 



2i\ R d 



24 dJ 2 k 



(i\ R ) 3 a d 



12 

(i\ R ) 3 a 



9Jk 



8 {Jl + 26)V2 



d 



{i\ R fa 



{Jl + 26)V2 



/ 



V V J * + 26 / 



(93) 



using ( 1621) and ( 1661) . This asymptotically behaves like 



ifcfc) fl =i ~ p 



(94) 



for k — > oo. Looking at (1331) . it is obvious that the higher-genus contributions are 
also decaying at least as rapidly as the genus 1 contribution (19^1) . and similar results 
could be derived for the n-point functions. 

Finally, it follows from (1631 that each term in F g contains the factor 



F g ~ ((i\ R ya 



(95) 



for g > 2, which is at least a 2 ^ s Hence {i\ R ) 2 a is the parameter which controls 
the genus expansion. 



5 Further aspects 
Structural remarks 

We briefly add some heuristic remarks which might shed new light on the formal 
results obtained above. 
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The model considered here is characterized by a function F of J which is invariant 
under conjugation with U(N 3 ). This implies in particular that F = F^(Jk) is a 
totally symmetric function of the eigenvalues Jk, i.e. it is a function on the quotient 
space M. n3 /Sn 3 where Sj^3 denotes the permutation group. Since the observables 
are obtained by taking partial derivatives w.r.t. J, we are interested in 

W k {J) = ^F N (J), 
oJk oJi 

etc. Renormalizability requires that all these derivatives exist and have a well- 
defined limit iV — > oo, hence that F/v(Ji) converges to an infinitely differentiable 
function F(Jj). This must hold at the point J with eigenvalues given by (fIBi) . 
or preferably in some "good" neighborhood U 3 J in R , with an appropriate 
limit N — > oo. Furthermore, these partial derivatives should have suitable decay 
properties such as indicated above. 

Now consider the dependence of Wk(J) on the index k. It is a physical re- 
quirement that this dependence on the indices is mild§. This can be understood 
analytically as follows. Note that Wk(J) is related to Wi(J) by exchanging the ar- 
guments Jk and Ji, i.e. by applying the permutation operator a^i on the space R^ 3 
of J k , 

W k = W l oa k ,i. (97) 

This means that the index dependence may be traded for a "small" change (a permu- 
tation) of the eigenvalues. Since the eigenvalues are ordered and approach a simple 
distribution, it is plausible that a^z respects the neighborhood IA. Smoothness in U 
would then naturally imply that the dependence on the indices is mild. Indeed, we 
found explicitly that the index dependence at genus becomes translated into the 
dependence of a smooth function W(y) f j43l on the variable y. 



Extension to Q ^ 1 

Once a suitable domain U is established where F[J] is smooth with suitable decay 
properties of its partial derivatives, then the following strategy to extend our results 
to Q ^ 1 can be envisaged. The partition function for ^ fl ^ 1 can be obtained 
from the results for Q — 1 using 

Z[J] = (e&rfcMZi,*])^ = e eTr[ ^ ][ ^ ] Z[J} n=1 

= e^i.^Pi,^] e F[J]n=i = . e F[J] (9g) 
8 some mild singularities for coinciding indices might be allowed however 
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for small e, accompanied by a change of wavefunction normalization. Hence we 
should consider the operator 

A 7 := Tr[% ^p,, -1] (99) 
oJ 2 oJ 2 

acting on functions of J. It does not commute with U(N 3 ), hence Z[J] will be a 
function of J which no longer depends on the eigenvalues only. We have established 
in this paper that Z[J}^ = i is a well-defined, infinitely differentiable function (after 
subtracting a possible infinite constant from F). The question of renormalizability 
for SI 7^ 1 is whether F[J] is also a well-defined, infinitely differentiable function, 
possibly after further renormalization. To establish this using the methods of the 
present paper might therefore be feasible, by studying the operator Aj and estab- 
lishing careful estimates on the partial derivatives of F[J]q = i and their asymptotic 
behavior for large N. However, we will not attempt to do this in the present paper. 



Remarks on the relation with string theory 

The Kontsevich model has been related to string theory in [19] as follows. The 
matrix Xij of the Kontsevich model is interpreted as coefficient of the open string 
field, more precisely the tachyon, connecting the (Liouville) D-brane with label i 
to the Z)-brane with label j. The eigenvalues Mi of the external potential in the 
Kontsevich model are interpreted as boundary cosmological constants on the brane 
with label i, and the Kontsevich model (1211 describes the (topological sector of) 
open string field theory in this situation. 

Applying this interpretation to our model, we are led to a picture of open string 
tachyons in string field theory with iV 3 D-branes, with specific cosmological constant 
J% ~ H + H + ^3 on the brane i = (ii, i^-, i^)- The fact that our model is renormal- 
izable implies that the correlators become essentially smooth functions of these 3 
coordinates ii,i 2 ,i3, with a well-defined large N limit. This could be interpreted 
by saying that the endpoints i,j of the open strings effectively live in 3 dimensions 
with coordinates i\, t2,h > 0, forming a 3-dimensional wedge with a potential 
determined by Jj. In other words, 3 extra dimensions appear from the string point 
of view by some kind of stringy "deconstruction" of dimensions. On the other hand, 
we have the interpretation as 6-dimensional non- commutative field theory. This 
points to interesting directions for further studies. 



6 Perturbative computations 

We write the bare action (jBJ) as 

S = Tr(l(J0 2 + 2 J) + ^0 3 -(zA)A0) 
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= TrfyiiGn^tf + ^^+SS. (100) 
The finite kinetic term {Gr) 3 ^ = \5\5^{Ji + Jj) defines the renormalized propagator 

a;.? = = sffiJ^ = s}6t l,{A f 2 } QV (101) 

hl 3 3 J l + J J 3 1 + 1+ (fJ? R + 3) 

corresponding to the finite (renormalized) matrix 

J\ ni ,n 2 ,n 3 ) = l6n d 9 2 (ji+ 3+ f R9 )\n 1 ,n 2 ,n 3 ) 1 (102) 

using the notation n = rti + n 2 + n 3 ([H]). Here Xr is the renormalized coupling con- 
stant, valid at a "scale" given by m, i.e. assuming that the indices of the interaction 
vertex approximately satisfy 

imjmk = m. (103) 

The counterterms are collected in 

5S = Tr ( - [i\ R )Z x A<f> + 1 -{5J<P 2 + 2 <5J) + ^ 3 ), (104) 

where 

5J\n 1 ,n 2 ,n 3 ) = I6n 3 6 2 [{Z— — ) + (Z - 1) J R j |ni, n 2 , n 3 ), 
= -Z^ A, 

= X R (Z X -1) (105) 

is part of the counter-term, and <5/i 2 = {^i 2 — ^i 2 R ). It is then easy to see that the usual 
power counting rules apply (with suitable extensions to the case of higher genus [20]), 
where N plays the role of the cutoff A 2 . This can be seen in the following 1-loop 
examples: 

3-point function In 6 dimensions, the interaction vertex requires renormalization, 
induced by the planar one-loop 1PI graph (without external legs) in figure [H This 
gives 

(MirikdiPi = + (-^) 3 E jn I jr jr I jr jn I jr ~ % ~Y 

i\r , 1 iSX 



2 +(-^zZ {J R + J R r - : Y+f Me 
i\ R (—iX R ) 3 sr^ 1 iSX 



(-tX R y ^ 1 idX 
2 + (16^)3 Q + 2 + 1 

iX R {-iX R f 1 N + m idX 
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Figure 1: one-loop contribution to 1PI vertex 



All nonplanar contributions are finite. This gives the counterterm to the coupling 
constant 

{\bir i U' L y 1 8m 



so that 

The minus sign indicates asymptotic freedom. However, the last formula (11081) is 
only suggestive, and such a one-loop result should be used with caution. Note that 
we could equally well write 

(<A) = - ^ - 2 + 0(X R ), (109) 

^ ^ 2 (16tt36I 2 )3 1u & m J 

which is in fact agrees much better with the exact scaling (1681) of the bare coupling. 
A better way to understand the coupling constant renormalization is to determine 
its dependence on the "scale" m, i.e. the running coupling constant 

g (m) = -i\ R (m) = 2{(j) ij (j) jk (j) ki )ip I (110) 

for i j « k = m (11031) . An extra factor 2 is inserted for convenience. The 
corresponding one-loop beta function can be obtained from (I106p . 

a - = 2m d ^jk4>ki)iPi = 9{mf d\og^ = g{mf 

P ' dm " (16tt 3 # 2 ) 3 Dm (16vr 3 # 2 ) 3 1 ' 

indicating asymptotic freedom. This gives 

dg(m) 1 dm 

g(m) 3 (167r 3 # 2 ) 3 m 

which can be integrated to give the running (1PI) coupling constant 



(112) 



g{mf = . (113) 

26 



This decreases with increasing scale m, which means asymptotic freedom. How- 
ever, note that the exact scale dependence g(m) ~ ~ji^mf f° r l ar g e m which was 
determined in (l§2"j) is not correctly reproduced by this one-loop RG computation. 
This shows that the common practice of using the 1-loop RG results for the running 
coupling constant may not be sufficiently accurate for large scales. 

It is worth pointing out that the coupling constant runs here even at one-loop, 
in contrast to the case of the 4 model in 4 dimensions [16]. 



1-point function The one-loop contribution to the 1-point function gives 

iX R Z x A i\ R 1 2 

A; 



hi) 1-loop ~ Ai j J J 



i\ R l „ . 1 1 



-Z x Ai + 



T 3()2 i h 4- 3 4- >i 2 -f) ' ( ' 114 ' ) 



J t \ 16vr 3 fl 2 + k + 3 + 

To proceed, we expand 

Hi) := J2 i + k + 3 + = m + (i) h ' (0) + 1 { - )2 h " (0) + hR ® (115) 

where hfj(i) is a finite nontrivial function of i, and 

"(o) = ErT^ = ( - 6l0g2+ 5 l0g3)Ar! 

+(-6(n 2 R + 3) log 2 + 3(/40 + 3) log 3) AT + -(/i 2 ^ + 3) 2 log N, 
h'(0) = - E tTT ^ T oN 2 = -(16vr 3 g 2 ) 2 / - (16vr 3 g 2 ) 3 (3 + 

= -(61og(2)-31og(3))iV + ^+^log(iV), (116) 
*'(0) = 2 £ ( , + ^ + 3)3 - -(16. 3 ^) 3 2/l = log(iV) (117) 

up to finite corrections, using the results of section 14.1.11 We note in particular 
that the i- dependent term (£) h'(0) + \(i) 2 h"(0) in ( 1114p forces us to introduce a 
corresponding counterterm to the action as in (jHJ), ^4 = a + cii J +02 J 2 . As discussed 
in section [21 this is no longer equivalent to an infinite shift fllOp of <fi. Taking this 
into account we have 

iX 



*>a)i-ioop = ~—f- ( - (oo + % J + a 2 J 2 )Z x 



+ w^ im + (i) ^ (0) + 1 & h " {0) + ' (118) 
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Finiteness and the condition (0oo) = implies to lowest order 

a, = ^^ys 1 2^N = -f 1 + fm i te, (119) 

ai = -(lQn 3 6 2 )(fi 2 R e + 3)a 2 + ^ 3Q2)2 h'(0) = -f + finite, (120) 

a ° = T77^^(0)-(87r 3 ^)(^ + 3)a 1 -(87r 3 ^) 2 (^ + 3) 2 a 2 (121) 

up to finite corrections. These renormalization conditions guarantee that the one- 
point function ((pa) has a well-defined and nontrivial limit N —>■ oo. In particular, 
PSD and dHJ) imply 



-= 1 / i +(^ i ^= i 4(^ 1o « jv ' < 122 ' 

to lowest order, and ( I120p together with (ITTj) gives 

^ = 5^J /o + /<n<te - (123) 

This is consistent with the exact result (1721) . 

Finally, it is interesting to consider the behavior of {(pa) for large i. After im- 
posing the above renormalization conditions, we observe that the dominating term 
in (1118ft is l -^rh{i) 2 h"(0). Here h"(0) ~ log(iV) should more properly be replaced 
by \og(N/i), which implies 

(0«>~ilog(i) (124) 
for large z. This is in agreement with (1821) . 



2-point function Next we compute the leading contribution to the 2-point func- 
tion (4>u4>kk) for Z ^ k, which vanishes at tree level. The leading contribution comes 
from the nonplanar graph in figure [21 which gives 




Figure 2: one-loop contribution to {4>u^>kk) 
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^ jkJi \Jk + Jij 

(for / ^ k) indicating the symmetry factors, where the disconnected contributions are 
given by (11141) . This is consistent with the result (1861) obtained from the Kontsevich 
model approach. Note that the counterterm 6 J does not enter here. 

Similarly, the leading contribution to the 2-point function (<fiki<fiik) for I ^ k has 
the contribution indicated in figure El involving also the counterterm 5 J. This gives 




Figure 3: one-loop contribution to 



>ki<Pik 



■>kl<Plk/l-loop 



J? + J? 

4 



2(iX 



R 



+ 



E 



(4>kk + 4>u) 



l 



SJi + 5J k 



(126) 



The first term is the free propagator, the next term the tadpole contributions, and 
the last them the one-loop contribution in figure [3] with counterterm 8 J. 

We have to adjust the parameters such that the result is well-defined and non- 
trivial. Using (11161) and (11171) . we can write 



E 



E 



k + l 



l 6vr 3^2)2 + fj^B + 3) 2 (j + fi 2 R 9 + 3Y 



+ finite 
(127) 



Therefore 



SJ, + SJ k 



(16tt 3 # 2 ) 2 
which using (11161) . (11171) implies 

(Z-l)(k + l) = 2(i\ R ) 2 (k + l)f 1 



(-/i'(0)-(fc + Z)~/i"(0)) + finite, (128) 



(129) 
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and 



Z(8/i 2 e) + (Z - 1)(3 + gf) = ^SLh + 2(^) 2 (3 + fifl) h (130) 



;16tt 3 # 2 )' 



up to finite corrections. Hence we obtain the lowest order mass and wavefunction 
renormalization: 

Z = l + 2 (zA i? ) 2 /i = - -A 72 



1 + 2 (16^02)3 10 8 JV 



2(iA^ 2 



= <iW° (131) 

up to finite corrections, in agreement with ( 112311 . 

Note that in 8 or higher dimensions, divergent mixed terms (k) (l_) would occur in 
( 112711 which can no longer be absorbed. Then the model is no longer renormalizable, 
as in the commutative case. 



7 Summary and discussion 

We have shown that the selfdual NC </> 3 model in 6 dimensions can be renormalized 
and essentially solved in terms of a genus expansion, by using the Kontsevich model. 
This provides a model which contains essentially the full complexity of renormaliza- 
tion of a not super-renormalizable asymptotically free quantum field theory, while 
being solvable and hence fully under control. In principle, all n-point functions can 
be computed in a genus expansion, and we give explicit expressions for the 1-, 2-, 
and some 3-point functions. 

In particular, we were able to determine exactly the RG flow of the bare pa- 
rameters as a function of the cutoff N, as well as the running of the "physical" 
coupling constant i.e. the 1PI 3-point function. As in the case of 2 and 4 dimen- 
sions [1,2], it turns out that the renormalization is fully determined by the genus 
sector. In particular, we can compare the exact results with the standard per- 
turbative methods. For example, it turns out that the 1-loop beta function for the 
coupling constant gives roughly the correct behavior and correctly predicts asymp- 
totic freedom, but wrongly gives a (log TV) -1 dependence on the scale as opposed to 
the correct (logiV) -2 dependence. 

We also show that the model has a critical surface defined by a = 0, which 
separates 2 different phases. One phase has the expected "physical" properties, 
while in the other some modes become unstable. 

It is very remarkable that a nontrivial, asymptotically free 6-dimensional NC <p 3 
field theory allows such a detailed analytical description. There is no commutative 
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analog where this has been achieved to our knowledge. Therefore this model can 
serve as a testing ground for various ideas and methods for renormalization. It also 
shows that the noncommutative world in some cases is more accessible to analyti- 
cal methods than the commutative case. While the techniques used in this paper 
are more-or-less restricted to the 3 interaction, it is worth pointing out that the 
renormalization is determined by the genus contribution only, which is accessible 
in a wider class of models; see also [13, 14] in this context. 

One of the open problems is the lack of control over the sum over all genera g. 
We have not shown that the sum over g converges in a suitable sense, which would 
amount to a full construction of the model. Furthermore, it would be extremely 
interesting to extend the analysis beyond the case of Q = 1. We propose a strategy 
in section [5] how such an extension might be possible, and we also comment on a 
relation with open string field theory. 
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